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Abstract 

We present a cellular automaton simulating the behavior of public bus 
transport in several Mexican cities. The headway statistics obtained from 
the model is compared to the measured time intervals between subsequent 
bus arrivals to a given bus stop and to a spacing distribution resulting from 
a random matrix theory. 
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The public transport in Mexico is organized differently to that known in Europe. First, no leading com- 
panies are responsible for the city transport. Thus, there are no timetables for the city buses and sometimes 
even not well defined bus stops. Moreover, the driver is usually the owner of the bus and so his aim is to 
maximize the income. Since every passenger entering the bus has to pay, the driver tries to collect the largest 
possible number of passengers. When not regulated the time interval under which two subsequent buses pass 
a given point will display a Poissonian distribution. This is a consequence of the absence of correlations 
between the motion of different buses. Such situation is, however, not welcomed by the drivers since then 
the probability density that two buses arrive to a bus stop within short time interval is large. In this case 
the first bus collects all waiting passengers and the second one that arrives bit later will find the stop prac- 
tically empty. This simple reasoning makes clear that existence of certain correlations between buses, that 
will change the Poisson process, will be of favor. Indeed, in Mexico various strategies have been developed 
to create such correlations. Here we discuss the situation in three cities: Cuernavaca, Puebla and Mexico City. 



In Puebla there is not an additive mechanism that helps to increase the bus correlation. Hence the head- 
way statistics of the buses in this city should be close to Poissonian. In Cuernavaca in turn, the information 
about the times, when the buses pass certain points, are notified and then sold to the bus drivers (there is a 
real market with this information) . The driver can in such a way change the velocity of the bus in dependence 
on the position of the bus in front of him and a bus-bus correlation appears. 
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To describe the correlations we use a modified version of the celebrated Nagel-Schreckenberg cellular 
automaton (see Jl|] and [||). Consider TV equal cells on the line (for our purpose we define the length of one 
cell as 30 meters) and n particles (buses) moving along it. Thus, c = j£ is the bus density and d — c^ 1 is 
the mean distance between two neighboring buses. Furthermore, define the maximal velocity v max of the bus 
and its probability p to slow down p G [0, 1]. The Nagel-Schreckenberg model describes the dynamics of the 
system with the help of the following update rules. In time t — the positions Xi of the particle i = 1 . . . n are 
integer numbers randomly chosen from the set 1 ... N satisfying the condition > Xi for every i = 2 . . . n. 
Furthermore, in time t = the velocity Vi of the z-th particle is set to zero for all i : Vi(0) = 0. The buses start 
to move with velocity v according to the update rules that has to be applied simultaneously to all particles. 

• Step one: When velocity of the bus is smaller than a maximal velocity 
Vmax, it increases its velocity by one 

Vi(t + 1) = min{^(t) + l;v max } (1) 

• Step two: Particles with positive velocities are randomly slowed down 

Vi(t + 1) =«i(t + l)-l (2) 

with probability p 

• Step three: Particles update their positions according to 

Xi(t + 1) = mm{xi(t) +Vi(t+ l);xi-i(t + l) - 1} , (3) 

i.e., the particles move according to the rule Xi(t + 1) = x.i(t) + Vi(t + 1) 
with the restriction that they cannot occupy the same cell or overtake 
each other. In that case the particle hops to the cell behind the occupied 
one. 

To create a modification of the model and to adapt it to the situation in Mexico we change the update 
rules on a subset M of possible bus positions where the information is passed to the driver. (The density of 
those points will be denoted by a; a = M/N.) 

To change the model we add to the steps (§]),(!]),(§) an additive step ([J]) that takes into account the 
processing of the information : 

At the points j G M the information about the time interval At to a preceding bus, which passed that 
point, is available to the driver. Using it the driver can modify the bus velocity - to speed up if At > t or to 
slows down for At < t. (Recall that £ is the mean time interval between subsequent buses.) Hence we change 
the model by adding 

• Step four: 

v j (t)=v j (t) + l 

if U(j)-ti-i(j) < 1 

v i {t)=v j {t)-l (4) 

if — ti-i(j) > i where U{j) denotes the time when a bus i passed 
through the point j: Xi{U{j)) = j. 
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Using the modified model we focus on the headway statistics, i.e. on the probability distribution of 
the time intervals At between two subsequent buses passing a given point and compare it with the results 
obtained in the cities. For the simulation of the transport in Puebla we choose c = 1/40, v max = 2,p — 0.5, 
and a = taking into account the fact that there are no checking points. It is not surprising that the headway 
statistics is in this case close to the Poisson distribution (see Figure 1) 

P(t) = e~\ (5) 

when t is the spacing re-scaled to the mean distance equal 1. Besides the headway statistics, we compare also 
a number variance T, 2 (t) that is defined as 

^(t) = ((n(t)-t) 2 }, (6) 

where n(t) is the number of bus arrivals to a given point during the time period of the length t. Note that 
(n(t)} = t due to the fact that (t) = 1. It can easily be checked that for a Poissonian process 

Y?{t)=t. (7) 
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Re-scaled time headway 



FIG.l. Time headway of the bus transport in Puebla. 
The curve represents the Poisson distribution (H). Plus signs 
display the time headway distribution of buses in Puebla and 
bars are taken from the CA model with a = 0. 

The number variance obtained from the data and from the simulation fits quite well with this prediction 
(see Figure 3). The data show, however, a small deviations from (Q). This is a manifestation of a weak 
interaction between the buses which probably originates from the fact that the buses interact through the 
number of passengers waiting on the stops. Namely when the distance between the buses is large, more 
passengers are waiting in the stop and the delay of bus in the stopping-place is longer. 

To simulate the situation in Cuernavaca and in Mexico City we use parameters c = 1/40, v max =2,p = 
0.5, a = 1/36, which represents one bus per 1.2 kilometers and one check point per one kilometer ap- 
proximately. The modified cellular automaton leads in this case to significant changes in the time interval 
distribution (see Figure 2) . The distribution obtained from the automaton fits well the observed time interval 
distribution. Moreover both distributions conform well with the distribution (see Ref. 
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P(t) = -^t 2 e^ 



(8) 



that describes the spacing distribution of a Gaussian unitary ensemble of random matrices (GUE). (It is 
known that this function describes also the distance distribution of certain one-dimensional interacting gases 
(see Ref. |)). 

Similar agreement is observed also for the number variance (see Figure 3) where GUE leads to 



with 7 re 0.57721566 (see (§). 



£ 2 (i) 



1 



(ln27rf + 7 + l) 



(9) 
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FIG. 2. Time headway of bus transport in Cuernavaca and Mex- 
ico City. 

The curve represents the Wigner formula (H.) Plus signs display 
the time headway distribution of buses in Cuernavaca and bars 
represent the results of the CA model with a = 1/36. The results 
obtained for buses in Mexico City are presented on the inset. 

However, as evident on the Figure 3, the interaction between the buses in Cuernavaca and in Mexico City 
is stronger than that resulting from the cellular model and leads to the stronger correlations. Whereas in the 
automaton model the correlations exist between the nearest neighbors only, in Cuernavaca and Mexico City 
one can recognize that interaction exists also between the first, second and third neighbor. 

We conclude that the modified Nagel-Schreckenberg cellular automaton successfully describes microscopic 
properties of the bus transport in some Mexican cities. The velocity of the buses is influenced by the informa- 
tion about their mutual positions so that the drivers can optimize their rank in competing on the passengers 
to be transported. This finally increase the coordination of the bus motion and changes the time headway 
statistics. 
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FIG. 3. Number variance of the public transport in some Mexican 
cities. 

Curve and line represent the prediction of the RMT and QVj) , 
respectively. Plus signs and stars display number variance ob- 
tained from the public bus transport in Cuernavaca and Mexico 
City respectively. Circles are the results of the CA model for 
a = 1/36. The number variance obtained from the transport in 
Puebla is displayed on the insert (plus signs). The circles show 
the result of the CA model obtained for a = 0. 
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